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Abstract 

The non-perturbative curvature inhomogeneities induced by relativistic viscous fluids are 
not conserved in the large-scale limit. However when the bulk viscosity is a function of the 
total energy density of the plasma (or of the trace of the extrinsic curvature) the relevant evo¬ 
lution equations develop a further symmetry preventing the non-linear growth of curvature 
perturbations. In this situation the fully inhomogeneous evolution can be solved to leading 
order in the gradient expansion. Over large-scales both the acceleration and the curvature 
inhomogeneities are determined by the bulk viscosity coefficients. Conversely the shear vis¬ 
cosity does not affect the evolution of the curvature and does not produce any acceleration. 
The curvature modes analyzed here do not depend on the choice of time hypersurfaces and 
are invariant for infinitesimal coordinate transformations in the perturbative regime. 
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1 Introduction 


The temperature and the polarization anisotropies of the cosmic microwave background are 
adequately described in the framework of the linear theory (see, e.g. the WMAP or WMAP9 
data essentially consistent with the Planck explorer data HI El El). However throughout the 
whole history of the plasma a non-linear growth of the curvature inhomogeneities cannot 
be excluded. The initial conditions of the temperature and polarization anisotropies of 
the cosmic microwave background are customarily set after neutrino decoupling pa ei hi 
but the evolution of curvature perturbations starts at a much earlier epoch namely during 
inflation and possibly even before. It is then particularly interesting, from the theoretical 
viewpoint, to define and discuss plausible non-perturbative generalisations of the curvature 
inhomogeneities in different physical situations going beyond the conventional cases of a 
single scalar field and of a perfect barotropic fluid. The theme of this paper concerns the 
non-perturbative generalization of the curvature inhomogeneities, their evolution and their 
physical relevance when the energy-momentum tensor is dominated by relativistic viscous 
effects at large-scales. 

Prior to the formulation of the inflationary paradigm the non-linear evolution of curvature 
perturbations has been analyzed in various frameworks. The approach followed here can be 
traced back to the expansion in spatial gradients of the geometry (see also 0). The 

gradient expansion has been used either in the proximity of the big-bang singularity or away 
from it. Close to the singularity the geometry may be highly anisotropic but it turns out to 
be rather homogeneous [SUB]. As soon as an inflationary event horizon is formed [HI El HU 
any finite portion of the event horizon gradually loses the memory of an initially imposed 
anisotropy or inhomogeneity so that the metric attains the observed regularity regardless of 
the initial boundary conditions as hypothesized in the past mmm- 

One of the central themes of the inflationary paradigm pans] is to wash out primeval 
anisotropies in the expansion right after the formation of the inflationary event horizon (see, 
however, Ref. mi). Probably the first non-linear generalization of inflationary curvature 
perturbations has been proposed in [9] after the pioneering analyses on the gauge-invariant 
treatment of linearised cosmological perturbations [10.J. These have been subsequently scru¬ 
tinized and rediscovered by different authors [18]. The non-perturbative generalizations 
discussed here may also have some impact on the neighbouring problems such as the higher- 
order approaches to cosmological perturbations (see |19] and references therein). 

It is appropriate to gauge the effects of the viscosity coefficients that can play a relevant 
role both in the early and in the late Universe. Since this investigation addresses the evolution 
of the non-linear curvature perturbation in the relativistic theory of viscous fluids (see, for 
instance, [20] and references therein) the assumption of the he strict reversibility of the 
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system will be dropped. The shear viscosity suppresses exponentially the traceless part of the 
extrinsic curvature. The bulk viscosity enters the definition of the curvature perturbations 
and may cause their non-conservations for typical length-scales larger than the Hubble radius. 
While the bulk viscosity does affect directly the non-linear deceleration parameter (possibly 
causing accelerated expansion) the opposite is true for the shear viscosity. 

From the technical viewpoint the variables introduced in the present analysis do not 
depend on the choice of time hypersurfaces and they are exactly invariant for infinitesimal 
coordinate transformations in the perturbative regime. It will be argued that the gravitating 
viscous fluids lead, under certain conditions, to a further symmetry preventing the non-linear 
growth of curvature inhomogeneities. The same set of conditions will also be shown to be 
compatible with a quasi-de Sitter stage of expansion. In the reversible limit our variables 
coincide with the ones conventionally defined in the case of a single scalar field or for a 
perfect barotropic fluid. The results pursued here can also be used to deduce the corrections 
to the linear theory without going through second or even higher-order calculations [ T9] , 

Finally, in the single inflaton case the large-scale cosmological perturbations are some¬ 
times treated within the so-called 5Af formalism [9] [18] where Af denotes the inhomogeneous 
expansion rate integrated in time generalizing to the non-perturbative level the total number 
of inflationary efolds. The presence of non-adiabatic fluctuations of the pressure make the 
formalism less appealing but we will show that, in some specific cases, the 5Af formalism 
can also be applied in the case of relativistic viscous fluid. 

The layout of the investigation is the following. In section [2] we shall discuss the ba¬ 
sic aspects of the geometry and of its coordinate transformations. The interplay between 
relativistic viscous fluids and general relativistic gradient expansion is addressed in section 
[3j Section [4] is devoted to the non-perturbative evolution of large-scale curvature inhomo¬ 
geneities in the viscous case. The perturbative limit of our considerations is discussed in 
section [5l Section E contains the concluding remarks. 

2 Nonlinear gauge transformations 

In the Arnowitt-Deser-Misner formalism [21] (ADM in what follows) the line element is 
expressed in terms of the conventional (3 + l)-dimensional decomposition: 

ds 2 = (? Mi /(r, x ) dx 11 dx v = N 2 dr 2 — ^{dx 1 + N l dr){dx :j + N^dr), (2.1) 

where N = N(t,x) denotes the lapse function, N l = N l (r,x) is the shift vector and 7 ^ = 
7 ij(r, x) is the spatial three metric E In the ADM variables of Eq. (12. 1 1) the extrinsic curvature 

2 Note that the spatial indices are lowered and raised using 7 ij(r,x). The Greek indices will take the 
values (/x, v) = 0, 12 3. 
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of the spatial slices (i.e. Kij(r,x )) and the components of the Ricci tensor of the spatial 
slices (i.e. r*j(r, x)) become: 


Kij(T,x) 


rij{r,x) 


1 

2 N 


—d-rlij + VjlVj + VjlVj , 


o (3)-pm «(3) pm _i_(3) pm (3)pn _(3) pm (3)pn 

^m 1 1 im ' 1 ij 1 mn 1 jn 1 im’ 


( 2 . 2 ) 

(2.3) 


where, for short, = V; is the covariant derivative defineqj with respect to the metric 7 ^; 
d T denotes a derivation with respect to the time coordinate r and ^ T™ are the Christoffel 
symbols computed from 7 ^. Note that T™' = ( T r”j but only in the case JV, = 0. From Eq. 
( 12 . ip the unit time-like vector normal to the x° = t = constant hypersurface is n^ = (JV, 0) 
and = (1/1V, —N l /N). The choices N — 1 and IV* = 0 corresponds to the geodesic slicing 
leading to the Gaussian normal coordinates, a popular gauge in numerical relativity [22], 
The condition IVj = 0 implies that coordinate observers coincide with normal observers: the 
normal vector has vanishing spatial component. Without positing a specific gauge choice, 
Eqs. (12.21) and (12.31) lead to the following expressions for the components of the Ricci tensor: 


Rq = 


R { = 


cLK 


N 


TrK + 


V 2 N N' 


N 


V m K H-. 

N N V 


VqK — VkKq ,, 


N 


V , I< - V k K[ 


1 1 ’ 

1 _, 1_ N m 


R] = —d r K] - KK:] - r] + - 


N T 1 '* ' N 

1 1 N j 

+ —V m W K™ - ViN m Kl - 

N N N 


N K ’ 

ViA'-ViA'f), 


(2.4) 

(2.5) 

( 2 . 6 ) 


where the notations K — 7 ' jl /V (( and TrA 2 — /\ ) A"' have been adopted. 

Notice that in the linear theory, for infinitesimal coordinate transformations of the type: 


are _> = x r + £ e > e e = ( £ 0 j £ <) j 


(2.7) 


the metric fluctuations change as the Lie derivative in the direction ra (see also section 
E3) • The perturbative expansion underlying Eq. (12.7j) assumes the separation of the geometry 
into a background value supplemented by a perturbation. The choice of the temporal gauge 
defines the spatial hypersurface of fixed coordinate time while the choice of the spatial gauge 
determines the worldlines of fixed spatial coordinates. The coordinate system is completely 
specified when both e° and e* are assigned. It is however possible to define various sets of 

3 We warn the reader that this identification will be followed throughout the paper. According to some, 
this notation may lead to potential ambiguities but we hope that, with this note, confusions will be avoided. 
In this paper V, will denote the covariant derivative on the spatial slices and not the spatial component of 
a (four-dimensional) covariant derivative. 
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gauge-invariant variables that do not change under Eq. (12.7j) . In the linearized treatment, 
the gauge parameters are of the same order of the metric perturbations: as soon as, in some 
gauge, the perturbation variables grow non-linear and affect the background geometry the 
linearised approximation is no longer tenable. 

The approach pursued here does not assume the validity of the perturbative expansion 
insofar as the geometry is not split into a background value supplemented by a perturbation 
with small amplitude: this treatment is arguably the most suitable for the unambiguous 
analysis of some backreaction problems. When the metric is not linearized around a specific 
background, the coordinate transformations will not necessarily be infinitesimal and shall be 
parametrized as = Y ,M {x) or, in more explicit terms, as: 

t->T = T(t,x), x* ->Y* = Y i (r,x). (2.8) 


Under the transformation of Eq. (12.81) the metric components of Eq. (12 .1 ft will change as 


can be written, schematically, a: 


4. 


The explicit form of Eq. 

(N 2 -N k N k ) = {a 2 -M k )(d T T) 2 -2p i d T Td T Y i -j ij d T Y%Y j , 


' IJ 

N, = -(a 2 - P t /3 k )d r TdiT + 2/3 k d r T d,Y k + 7 lk d T Y'd i Y k , 
Hi = -(c?-faP k )d j Td i T + 2p k d i Td j Y k +i lk d t Y k d j Y‘, 


(2.9) 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


where the lapse function, the shift vectors and the spatial three metric in the transformed 
system have been denoted, respectively, by a — a(T,Y), /3i — /3i(T, Y) and 7 ek = 7 ^( 7 , Y); 
Eq. (12. ip becomes, in the transformed frame, 

ds 2 = G^(T, Y) dY 11 dY v = o?dT 2 - 7 ^(dY* + pdT)(dY j + (3 j dT ). (2.13) 

We shall often refer to the concept of non-linear gauge-invariant variables. This terminology 
refers to the possibility of finding specific quantities that do not depend on the choice of time 
hypersurfaces and that are exactly invariant for infinitesimal coordinate transformations in 
the perturbative regime. 


3 Gravitating viscous fluids and gradient expansion 

Whenever dissipative effects are included both in the energy-momentum tensor and in the 
particle current the physical meaning of he four-velocity u M must be specified. In the Eckart 

4 The partial derivatives with respect to r and x l shall be denoted, respectively, by d T and by di. 
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approach u M coincides with the velocity of particle transport. Conversely, in the Landau 
approach the velocity u M coincides with the velocity of the energy transport defined by 
the (Of) component of the energy-momentum tensor giving the energy flux. The Landau 
approach shall be privileged mainly for practical reasons. 

3.1 First and second viscosity in the Landau frame 

The total energy-momentum tensor of the problem is given as the perfect held contribution 
(characterized by a total pressure p t and a total energy density p t ) supplemented by the 
irreversible contribution: 


T Z = C Pt + Pt)upU v - p t 6" + 7J( 77 ,0, (3-1) 

where denotes the viscous energy momentum tensor depending on the first and 

second viscosities (i.e. 77 and £): 

= 2<v = ±77 Vi W 7i . (3.2) 

As usual, Vp = (5p — upu v ) and the tensor Wyt appearing in Eq. (j3.2[i is defined as: 

WyS = VyUs + V gUy - - Pyg V\U X . (3.3) 

The total particle current will be denoted by = ntW M + zA (where n* is the total concentra¬ 
tion of the fluid while z/' denotes the diffusion current). Note that will denote hereunder 
the relativistic thermal conduction four-vector. As we shall see explicitly = 0 since 

can be written as zy = f(T,p t ,pt)V°d a JI where ]I = p/T and p is the chemical potential. 

The viscous energy-momentum tensor at large-scales can be evaluated in the Landau- 
Lifshitz or in the Eckart frames. In the Eckart case the four-velocity u M appearing in Eq. (13.11) 
denotes the velocity of the particle transport. The tota@ particle current vanishes in the 
comoving frame. Consequently the Eckart frame is fixed by requiring that j a u a = 0 while 
T^ u u u 7 ^ 0. The Eckart approach [23] (see also second paper of Ref. [2D]), seems preferable 
when the concentration of radiation quanta exceeds the concentration of the other species. 

Conversely, in the Landau-Lifshitz approach [23] pure thermal conduction corresponds 
to an energy flux without particles: the four-velocity u M coincides with the velocity of the 

5 In this paper we discuss the global evolution of large-scale curvature perturbations. We shall therefore 
deal with global quantities such as the total pressure, the total energy density of the system, the total 
viscosity of the fluid and so on. Of course the same analysis can be extended to the case where the fluid is 
composed by a number of fluids interacting among themselves as it happens, for instance, prior to photon 
decoupling. 
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energy transport implying T^ v u v = 0 . The two approaches are largely equivalent but the 
Landau-Lifshitz approach seems slightly more convenient, in the present context. Both the 
Eckart and the Landau approaches are suitable for the present class of problems where the 
typical scales are much larger than the mean free path. The second-order dissipative effects 
[23] become particularly relevant in the collisions of heavy ions |26j where, however, not all 
the numerous second-order terms have been so far included in the available analytical and 
numerical discussions. 

The covariant conservation of the total energy-momentum tensor (i.e. V^T 7 ^ = 0) can 
be projected along u v and along P“; the two obtained equations together with the covariant 
conservation of the particle current are given hereunder: 

+ (h)u*\ - u a d a p t + upVcT* = 0, (3.4) 

(pt + Pt)u 0 Vpu a - d a p t + u a upd 0 p t + P“V #1 T' M ' = 0, (3.5) 

V a {n t u a + u a ) = 0. (3.6) 


Using then Eqs. (13.41) and (13.61) together with the first principle of thermodynamics, the 
evolution of the entropj@ can be easily derived: 

X7 a [su a - pv a ] + u a d a JI = X7 a up T a0 /T. (3.7) 


In Eq. (13.71) . as already mentioned after Eq. (13.31) . p = p/T is the chemical potential rescaled 
through the temperature, s is the entropy density and u a is given by: 


= X 


n t T 

Pt+Pt 


d a /i — u a u 0 dpp 


(3.8) 


where x denotes the heat transfer coefficient. From the definition of the viscous energy- 
momentum tensor of Eqs. (13.2j) and (13.31) we can also explicitly write the term at the right 
hand side of Eq. (13.71) 


(V aU/3 )T a0 /T= (£/T)(V a u a ) 2 + 2(ri/T)a liU a^. (3.9) 


The adiabatic limit is recovered when the viscous contributions are neglected and the total 
entropy four-vector is conserved. The right hand side of Eq. (13.91) is positive semi-definite 
provided £ and p are both positive semi-definite. 

It is appropriate to remark, at this point, that the perfect fluid contribution is charac¬ 
terized by the barotropic index w = pt/Pt and by the related sound speed c 2 st = p'Jp'. In 
linear theory the fluctuations of the total pressure of the fluid are customarily decomposed 

6 The explicit form of Eq. (13.711 has been obtained by trading the term u v S7 for (V m m i ,)T M 1 ' since, 
in the Landau frame, V M (M y T Ml/ ) = 0. 
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into an adiabatic component supplemented by the entropic (or simply non-adiabatic) con¬ 
tributions 5pt = c 2 st 5p t + Jpnad (see, for instance, [57] and the Eqs. (j5.5jMj5.6j) in section 
[5]). This occurrence would correspond, at the level of the non-linear discussion, to the case 
w 7 ^ c 2 st where w may be a space-time dependent function (see section HP . Even if the con¬ 
ventional terminology might suggest otherwise, the non-adiabatic modes have nothing to do 
with the global viscosity of the system and may even arise in a globally inviscid fluid. This 
potential confusion of the standard terminology should be borne in mind to avoid unwanted 
misunderstandings. 


3.2 Gradient expansion of the Einstein equations 


The Ricci tensor reported in Eqs. (12.4[) . (12.5p and (12. 6 p have been already written in a 
form where the terms containing spatial gradients are distinguished from all the other. The 
same criteria must be adopted when expressing the explicit components of the total energy- 
momentum tensor so that, at the very end, we shall be able to write down the Einstein 
equations in their contracted form: 


K = 


rpb 


T ' 
-8 V U 
2 M 


ip — Vo7 tG, 


(3.10) 


where T = T“ is the trace of the total energy-momentum tensor and must not be confused 
with the effective temperature of the fluid appearing in the previous subsection. 

The parameter counting the gradients can be formally indicated as the gradient itself in 
units of the trace of the extrinsic curvature, i.e. as A = V//l(r, x). From Eq. (13 .1 p and 
bearing in mind the right hand side of Eq. (I3.10K we have: 


rri P ( i T~) \ Pt P zff 

Tfiv — (Pt p f) U H U v 2 


+ V <Vfj,u u + VvUfj, - u° 


oP^v "h Ujy V 


- -V IW V\U 


(3.11) 


The trace T has been expressed in terms of the effective pressure P e ff defined, in our case, 
as: 


T = T a = Pt~3 p eff, Peff = Pt ~ a U° . (3.12) 

The term V a u“ can be easily expanded in spatial gradients and the result is: 

1 c) 1!^ 

V«u“ = -K - —d k [N^u k ] + + 0(A 3 ). (3.13) 


From Eqs. (j2.4(1 . (12. 5 p and (12. 6 p . the terms containing the shift vectors turn out to be 
0( A) (see also [9] [IB]). It seems therefore appropriate to select the gauge N t — 0 where 
the coordinate observers coincide with normal observers. Equivalent choices, for the present 

























purposes, include the coordinate system where the expansion is uniform (i.e. K = K(r), in 
our notations) or the gauge where the energy density is uniform. 

The full expression of the various components of Eq. (13. lip are necessarily lengthy so we 
shall just exemplify the ( 00 ) component and then mention the leading order results for the 
other components. The (00) component of Eq. (13.111) is given by: 


T N 2 

T 00 ~ 9oo — + 3Peff) + N 2 u 2 (p t + P e ff ) — 2r/Nu 2 J-(N, u k ), (3-14) 

where J-"(N, 7 ^, u k ) is defined as: 

u k dkV 1 + u 2 + u k u^Kkj 

-_[d T (^VlT^) - d k (N^u k )]). (3.15) 


P(N, 7 ij, u k ) = \ d T [Vl + u 2 } - Uid l N - NVl + u 2 
1 


3^7 l 


Except for the first term at the right hand side, all the remaining contributions are 0( A 2 ) 
since they contain, at least, two spatial gradients. From the lowest-order form of the mo¬ 
mentum constraint (see below Eqs. (13.16jl and (j3. 191) ) the leading contribution of the spatial 
part of the velocity is clearly 0( A) so that u 2 = ^pjU l u = 0( A 2 ). 

The same procedure outlined in the case of the (00) component of the energy-momentum 
tensor can be discussed for the remaining components and the leading order results are: 


T° = t P, + Peff)u,u 0 + 2 V u t u 0 K t i +O (\ 3 ), 
= - 1 -(p t -P eff )S!-2r,Ki + 0(X 1 ), 


(3.16) 

(3.17) 


where the traceless part of the extrinsic curvature K] = Kj — 5{K /3 has been explicitly 
introduced. 

Using the results of Eqs. (13.141) . (13.161) and (I3.17|) together with Eqs. (12.41) . (12.5ft and 
(12.61) . Eq. (13.101) will become: 


—<%K - Tr K 2 
N T 

V i K - Vfc/vf = 
4 drKt - KK{ 


'l 


- y(P* + 3 Peff), 


N£i 


(p t + P e ff)uiU° + 2rjK\ujU l 


-r\ = e p 


(Pg " 2 — S{ - 2r)K‘ + n| 


(3.18) 

(3.19) 

(3.20) 


where III denotes the anisotropic stress which is by definition a traceless rank-two tensor in 
three dimensions. We have chosen to keep generic the form of fl^ since it may contain all the 
potential sources of anisotropic stress not necessarily connected to the the fluid sector such 
as scalar fields or even gauge fields. In all cases n) contains at least two spatial gradients. 


9 




































The anisotropic stress and the curvature rj are of higher order in the gradients but have 
been kept for the benefit of the forthcoming discussion aimed at showing the the traceless 
part of the extrinsic curvature is of higher order in the gradients and it the only component 
affected by the presence of shear viscosity. 

Finally, the evolution equations stemming from the covariant conservation of the energy- 
momentum tensor and of the particle current of Eqs. (13.41) . (13.5ft and (I3.6jl are given by 


- K ( p t + P e ff ) — 0 


—d T u l + u k 


N 

1 

N 


d T P f 


r r eff 


l N( Pt + p eff y k 
d T n t — Kn t = 0 . 


SI - 2 Kl 


d l N &P eff 
N P+Peff 1 


(3.21) 

(3.22) 

(3.23) 


The absence of the dissipative coefficients arising in the diffusion current is justified since 
these terms are of higher order in the spatial gradients, as it can be easily appreciated from 
Eq. (US]). 


3.3 Decoupling of the shear contribution 

From Eq. (13.201) the shear contribution only affects the evolution of the traceless part of the 
extrinsic curvature K\ and does not enter the deceleration parameter whose sign is solely 
determined by the bulk viscosity coefficient. 

Indeed, after taking the the traceless part of Eq. (13.20}) the following equation is obtained: 

d T K\ - NKK\ = -2riNf P Kl + Nf P Ilj + Nrj (3.24) 


where r{ — rj — Sj r/3 is the traceless part of the spatial curvature. In the general situation 
where 7](r,x) Eq. (13.24j) implies 




,f) 


£) e -2A(T*,T,g) 




v^( r ^ 

,x) 






+ 

f P 

1 


r dr" 1 


x)N(r", 

x) e~ 2A P "’ T ’ 5) n] 

(r", x)dr' 


V^( r ’ 

x) 

Jr* 

+ 

1 

/ 


r dr", 


,x)N(t", 

x)e- 2A(T "’ T ’ g) rjl 

{r",x)dr" 


V7(f 

x) 

Jr* 


(3.25) 


where 7 = det( 7 ij). I 11 Eq. (I3.25P r* = r*(x) denotes some arbitrary integration time while 
for two generic times 77 and r 2 , *4 .(ti, r 2 , T) is defined as: 


«4 (ti, t 2 , x) = f P [ 77(r', x) N{t', x) dr'. 

Jri 


(3.26) 
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Equations (13.25j) and (I3.26p show that the traceless part of the extrinsic curvature is deter¬ 
mined by the anisotropic stress and by the traceless part of the intrinsic curvature. Both 
quantities are of higher order in the gradient expansion^. Equation (13.26ft shows that the 
shear viscosity suppresses the traceless part of the extrinsic curvature even further in com¬ 
parison with the case 77 —> 0. The features of the damping are determined by the explicit 
expression of 77 . For a system dominated by radiation 77 ~ £ m f P Pt where £ m f p denotes the 
mean free path (for instance the Thomson mean free path prior to photon decoupling). In 
this case A — K£ m f p <C 1. In more general terms, however, 77 can depend on p t on the trace 
of the extrinsic curvature, on the total particle concentration and in all these cases Kij may 
even be much smaller than the anisotropic stress. 

Since 77 decouples from the trace of the extrinsic curvature, it does not contribute to the 
inhomogeneous generalization of the deceleration parameter. For the sake of comparison 
with the fully homogenous case we choose Gaussian normal coordinates and set N — 1; in 
this situation Eq. (I3.18p can be written as: 


q(t , x)TtK 2 = £] 


(Pt + Peff)UoU + 


0 , P eff ~ Pt 


27/7i 2 J'(l,7 ii ,M fc ), 


(3.27) 


where q(x, t) = — 1 + K/TyK 2 is the inhomogeneous generalisation of the deceleration pa¬ 
rameter] and the overdot denotes the derivative with respect to the cosmic time coordinate t 
which coincides with r in the case N — 1. The function .F(l, 7 ^, it fc ) (defined in Eq. (j3.15D ) 
accounts for the higher-order corrections. In general terms TrA ' 2 > K 2 /3 > 0, where the 
sign of equality (in the first relation) is reached in the isotropic limit. Since y *- 7 8 is always 
positive semi-definite, it is also clear that u^u 0 = 1 + Y^UiUj > 1. From Eq. (I3.27P it also 
follows that q(t,x) is always positive semi-definite as long as (p + 3 P e ff) > 0. This means 
that the sign of the generalized deceleration parameter only depends on P e ff (and hence on 
the bulk viscosity) while the shear viscosity does not play any role. According to Eq. (13. 2 7ft 
the correction of the bulk viscosity only arises to second order in the gradient expansion 
where, however, the bulk viscosity also contributes through the term (p + P e ff)u 2 implicitly 
contained in (p t + P e ff)uoit°. 

7 The higher order terms in the gradient expansion can be computed by following iterative methods where 
the spatial geometry is reconstructed, order by order, starting from a seed metric that do not contain any 
spatial gradient mm but this is not our primary goal in this investigation. 

8 In the homogeneous and isotropic limit, 7 ^ = a 2 (t)<5,j, K{ = —HSj and, as expected, q(t) —> —aa/a 2 . 
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4 Gauge invariant variables and their evolution 

4.1 Generalities 

Under Eq. (12.81) the energy-momentum tensor transforms as 

T^t,x) ->7V(r,?) = (^)(^)r^(r,x), (4.1) 

where T is the transformed time coordinate and will not be confused with the trace of the 
energy-momentum tensor. Equations (12. lOj) (j2. 1 1[) and (j2.12[) hold for the metric; similar 
expressions hold for the transformed components of the energy-momentum tensor. The 
analog of Eq. (I2.1(jp will be reported to fix the notations 

T 00 (t, x) = (d r T) 2 T 00 (T, Y) + 2(d T T)(d T Y i )T 0l (T, Y) + (<9 T V*)(d T r^(T, Y). (4.2) 

The explicit expressions of the (Oi) and (ij) components can be easily written in terms of 
the notations of Eq. (j4.2j) and will be employed below. 

The coordinate transformation must preserve the order of the gradient expansion. This 
implies, from Eq. (12. lip , that /3, = N, t — 0 and the coordinate transformation demands: 

a 2 d T Td{F = ^f ek d T Y e d 7 Y k . (4.3) 

The transformations preserving the order of the gradient expansion f5j can be written as 
follows: 

t^T = T(t, x), x k -> Y k (r, x) = / fc (r, x) + F k (r, x). (4.4) 

By construction the function f l (x,r) does not contain any gradient while F l (r,x) contains 
at least one spatial gradient; f k (r, x) can then be parametrized as f k (r, x) = x k g(r, r ) where 
r = \/xiX l . Since in the transformation all the spatial gradients of g(r, r ) will automatically 
contribute to F k {r , x), the effect of g(r, r ) boils down to a redefinition of a in the transformed 
frame. For this reason and for the sake of simplicity we shall set g(r,r) = 1. 

Equations (I2.10p and (I2.12p . thanks to Eq. (14 .3 p . will then lead, respectively, to the 
following pair of conditions: 

N 2 = a 2 (d r T) 2 - 7 ^drY'drYi, (4.5) 

Hi = -cPdiTdjT + T k ediY k djY e . (4.6) 

Recalling the explicit form of Eq. (14.40 . to lowest order in the spatial gradients, Eq. (14.6p 
implies that 7 ij(r,x) = 7 ij(T,Y) while Eqs. (14. 3 p and (14. 5 p determine the explicit form of 
F k (T,Y ); the explicit results are: 

F k (T,Y) = j <ffiV 2 (T)--^J, 7 „(r,x) = 7 «(T,F). (4.7) 


12 



























Equations (11.11) and (I4.7|) can be inserted into the various components of Eq. (14.ip to obtain 
the transformation properties of the pressure, of the energy density and of the velocity: 


P(r,x) = p(T, Y), P ef f(r,x) = P eff (T,Y ), 
Ui = Hi + adiT, N = a d r T. 


(4.8) 


4.2 Non-linear curvature inhomogeneities 


In the viscous case the non-linear generalization of the curvature perturbations on comoving 
orthogonal hypersurfaces is: 

Kx{t,x) = ^ Vi [hi (V7)] - 7^<9 r [ln (V7)] u i- ( 4 - 9 ) 

Using the transformation properties defined by Eqs. (14. 3 j) . (14. 5 j) . (14. 6 j) and (14. 8 j) we have 
that transforms as: 


Ui(T,x) -)• 


v (T V\ - 

’ J 3 dYJ dx i 


1 (_ 

—— I Hi “f" ct 


3a 


dT \ cl[ln (yq)] _ 
dxw dT 


1 ^[ln (-y/r)] dT 
3 dT dx i 


(4.10) 


since the two intermediate terms simplify in Eq. (I4.10p . we have that the curvature inhomo¬ 
geneities are invariant i.e. lZi(T,x) = lZi(T,Y). 

Using the same strategy applied in the case of Eq. (14.91) . the non-linear generalization of 
the density contrast on uniform curvature hypersurfaces becomej^: 

C,(r,f) = lv i [ln(^y)] + 3(p ^ ) . (4.11) 

The same analysis leading to Eq. (I4.10p can be performed in the case of the (i(r,x): 


C 


->■ 

+ 


7 (T y , = 1 d i ln (Vt)] d Y 3 1 ^[ln (Vt)] dT 

; 3 dY'J dx i 3 dT dx i 

1 dp dY* 1 dp dT 

3 (p t + Peff ) dYi dx i + 3(p t + P ef f) dT dx i ' 


(4.12) 


Equation (13.2ip (stemming from the covariant conservation of the total energy-momentum 
tensor in the transformed frame) implies that the derivative of ~p t with respect to T equals 
aK(j) t + P e fj). Thus the second and fourth terms at the right hand side of Eq. (14.12p cancel 

9 In linear theory the density contrast on uniform curvature hypersurfaces is invariant under infinitesimal 
coordinate transformations. Since, by definition, it has the same value in different gauges it can be also 
interpreted as the curvature perturbation on the hypersurfaces where the energy density is unperturbed. 
These two physical interpretations are relevant when discussing the so-called SAf formalism (see section (S]). 
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since aK = —9 In y/^y/dT. As in the case of Eq. (14. lOj) the invariance of Q is manifest since 

Ur,x) = UT,Y). 

Since, in the general situation, the bulk viscosity coefficient is a space-time scalar function 
its derivative transforms non-trivially under Eq. (1 1. Hi : 

f)f d£ — df 

^ = w +dr(Kr ’ « = ^ r - (413) 

Thanks to Eq. (14. 13ft we can obtain a further non-linear variable invariant under Eqs. (14. 3 ft — 
(14.5 j) and (I4.6]) - (I4.8I) that has no analogue in the inviscid case: 

1 KN 8 £ 

Z,{t,x) = -Vjln(^)] + —(4.14) 

From the gauge-transformed expression of Eq. (14.141) . using Eqs. (14. 5 p and (14.61) the gauge- 
invariance of Zi(r,x ) is easily demonstrated. 

Let us finally mention, for the sake of comparison, that the general form of Eq. (14.10)1 
can be used to recover the well known results obtainable in the case of the single scalar field 
[9]. The (Of) component of the energy-momentum tensor of a minimally coupled scalar field 
= di(pd r (p/N 2 implies that tq = NVi<p/(d T (p). Therefore Eq. (14.91) implies 

Ki(T,x) = ^Vj[ln (V7)] + ( 4 - 15 ) 

which is also invariant since, from Eq. (14.8f) . the transformation of <p will be given by: 

d T ip = d T Jp(d T T ), -^r = + 8 T TpdiT. (4.16) 

Equations (14.15j) and (14. 16j) reproduce the standard results of Refs. [2, 18j. 


4.3 Evolution of the gauge-invariant variables 


According to the momentum constraint of Eq. (13.191) the combination Kut appearing in Eq. 
(14.91) is expressible in terms of the gradients of the extrinsic curvature as: 


Kui 


1 

3 P P (Pt + Pel,) 


Si I< 2 - 3A-VU4 


(4.17) 


In Eq. (14.171) and in the forthcoming discussion we shall keep the dependence on the traceless 
part of the extrinsic curvature just to keep track of the difference between and Q. From 
Eq. (13.18[) and from the trace of Eq. (I3.20P we obtain the following pair of equations: 


2 f P p t = K 2 - Tr K 2 , 3N£ 2 P (p t + P eff ) = 2d T K - 3 NTtK 2 + NK 2 , (4.18) 
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implying that Eq. (I4.17P can be finally expressed as: 


Kui 


dipt 

Pt + Peff 


+ Gil 


Gi = 


6 P P {pt + Peff) 


3d ?; [TrA 2 ] -GKVkK 1 ; 


(4.19) 


After inserting Eq. (14.19ft into Eq. (14.9ft . the partial time derivative of P % becomes: 

a T n : = -\dANK) + \dJ d ‘ p ‘ ) + (4.20) 

3 3 \pt + P e ffJ 3 

To leading order in the spatial gradients, the same kind of evolution equation reported in 
Eq. (I4.20p is derivable for Q starting directly from the definition Eq. (14.lip and using Eq. 
(14.18ft . The leading terms of the evolution equation will be the same and the rationale for 
this occurrence is that Q an d Pi differ by terms that are of higher order in the gradient 
expansion, i.e. 6 i 2 P (Q — Pi) = — [3<9j(TrA~) — 6 K'V^K^]. 

Inserting now Eq. (13.2ip into the first term at the right hand side of Eq. (14.20ft we arrive 
at the following result: 


d T Pi = -d 7 


/ dipt \ 

^ d 

( d T p t \ 

\Pt + Peff) 

3 d * 

\Pt + Peff) 


+ 


d r Gi 


(4.21) 


The third term at the right hand side of Eq. (14.21ft will now be dropped since it is of higher 
order in the gradients. Equation (14.21ft can be expressed in a physically more significant 
form by separating the viscous contributions from the conventional non-adiabatic terms 
that normally appear even in the absence of irreversible contributions: 


d T Pi 


S na d (A P) T ^viscous (A P)i 


(4.22) 


where the two source terms are given, respectively, by 


S na d ( T , x) 


KN 


-(dipt - c 2 st dip t ), 


3 (Pt + Peff) 

K 

&viscous , x') 0 / | \o\(drPt)di£ i diPt{d T S))\ 


3 (pt + Peffy 


(4.23) 


+ o/ Jp \2 [(d T Pt)diK — {dip t )d r K}. (4.24) 

3 [pt + PeffY 

As already mentioned, the total sound speed is c 2 t = d T p t /d T p t . Let us consider first S na d(T, x) 
and show that it is nothing but the standard adiabatic contribution. Broadly speaking the 
barotropic index is a space-time function, i.e. w(t,x ) = Pt/Pt and Eq. (I4.23P implies: 


S n ad (j x'j 


KN Pt 


3 (Pt + Peff ) 


di.w - 


PtdiPt 


3 (Pt + Peff ) 


-d T w. 


(4.25) 


In linear theory when w is a space-time constant the sound speed equals y/w and the non- 
adiabatic contribution is absent. In this situation S na d(r, x) —> 0, as expected. 


15 




































According to Eqs. (14.221) . (I4.23P and (11.2 111 the relativistic viscous fluids lead to a source 
term implying that 77; is not constant in general terms. This may happen in various situa¬ 
tions where, for instance, the viscosity coefficients depend in time and in space. For instance, 
across the matter-radiation transition the shear viscosity coefficient rj determines the optical 

depth, the Silk damping scale and, ultimately, the shape of the visibility function [4] [28]. 

/ \ 2 <? 

The scaling properties of q and £ can be expressed £/r? ~ | — c £ t ) for q > 1. Across the 

matter-radiation transition the sound speed of the plasma interpolates between l/\/3 and 0. 
Adopting the viewpoint of linear theory (see section [5]) and separating the background from 
the fluctuations the sound speed can be computed as c 2 t = 4/[3(4 + 3a)] where a = a/a eq 
denotes the scale factor normalized at equality. This dependence in £ implies the generation 
of non-adiabatic modes (see, in particular, the last paper of [27]). 

4.4 Conservation of curvature perturbations 

If c 2 t Y w (or if w is a space-time function) the curvature perturbations are non conserved 
even in the limit £ —» 0, and this is nothing but the standard situation of the conventional 
non-adiabatic modes. To exclude all the potential sources that could make 77,; time depen¬ 
dent besides the ones we ought to investigate specifically, namely the relativistic viscous 
contributions we shall posit that w is constant and that c s t = \/w. This choice implies, 
according to Eq. (14.231) . that S na d = 0. 

The only contribution remaining at the right hand side of Eq. (I4.22p is the one com¬ 
ing from S viscous . Furthermore, by focussing on Eq. (14.241) we see that the second term 
(proportional to £) can be rewritten as 

—i—j[(S Tft )5.A- - (d iPt )d T K] = (A d,(K 2 - 3{%p t ), (4.26) 

3 (Pt + PeffY 6 (p t + P eff ) 

but the term at the right hand side vanishes because of the first of Eq. (I4.18P : in fact, 
Tr K 2 = K 2 /9 + TrK and TrK is of higher order being proportional to the square of the 
total anisotropic stress. We stress that the result of Eq. (14.261) holds non-perturbatively; 
it does not assume a separation between the background space-time and its perturbative 
fluctuations. Equation (14.221) becomes then: 

drlli = o/ \ w [( d rPt)diZ ~ dMd T a (4.27) 

3 (Pt + Peff) 2 

If the source term in Eq. (14.271) vanishes the curvature inhomogeneities will be conserved 
and the equations of motion will enjoy a further symmetry since 1Zi{r,x) can be shifted by 
a a term constant in time (but not in space). The bulk viscosity £ can depend, in principle, 
on five quantities, namely K, TrK 2 , p tl p t and n t . The dependence on p t can be traded 
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for p t since p t = wp t with constant w. Since the traceless part of the extrinsic curvature is 
of higher order we can also drop the dependence on Tr K 2 that coincides, to leading order, 
with K 2 /9. Equation (I4.18P can be finally used to relate K 2 and p t . To lowest order in the 
gradient expansion we have therefore only two qualitatively different cases: £ = £(p t ) and 
£ = £(p t , n t ): the case £ = £(n t ) is indeed the same as the one where £(p f ,n*). 

When £ = £(p t ) Eq. (14.271) implies that d T lZi = 0: in this case the two terms at 
the right hand side simplify because d j£ = (<9£/ dp t )dip t and <9 T £ = ( d^/dp t )d T p t . Putting 
together the results obtained so far, we can therefore say that Eq. (14.211) is invariant for 
77.,;(r, x) —* 7 Zi(r,x) + Q(x) provided £ is either a space-time constant or a function of the 
total energy density. 

The requirements of the previous paragraph correspond to the situation where Eq. (14.181) 
admit a fully inhomogeneous solution whose homogeneous limit is of quasi-de Sitter type. 
Let us therefore show explicitly that this is indeed the case in the simplest situation where 
£ is a space-time constant. Equations (14.181) implies the following decoupled equation^ for 
I<: 

K - ^AT 2 = | tliK, (4.28) 

The general solution of Eq. (j4.28[) is: 


K(t, x) = 


K 0 (x) e K ^- to ^ 2 


(w + l)K 0 (x) 


X _ e K i [t-t. 0 {x)\/2 


+ K £ 


IU = 3 ^£ = 


V 


(4.29) 


where Kq(x) constant in time but not in space and is just a parameter of the solution. 
For K^[t — t 0 (x}} C 1 we have that K(t,x) is singular while in the opposite limit it goes to 
a negative constant. In the homogenous limit to is constant also in space. 

Equations (14.281) and (14.291) show that when £ is a space-time constant we can derive 
an inhomogeneous solution whose homogeneous limit interpolates between a perfect fluid 
solutions and a quasi-de Sitter solution. These solutions are the inhomogeneous counterpart 
of various quasi-de Sitter solutions derivable in the fully homogenous limit (29j. The gener¬ 
alization of Eq. (14.28P to the case when £ = £(p t ) is straightforward since the dependence 
on p t can be eliminated through Eq. (14.181) . Some cases where £ has a power-law depen¬ 
dence on p t are even analytically solvable. We can therefore conclude that the large-scale 
curvature inhomogeneities are non-perturbatively conserved in numerous cases where the 
non-perturbative solution, in its fully homogeneous limit, admits a set of quasi-de Sitter 
backgrounds. 

10 In Eq. (14.281) we have chosen the geodesic slicing with N = 1; in this case t = t where t denotes the 
cosmic time coordinate and the overdot denotes a derivation with respect to t. 
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4.5 Non-conservation of curvature perturbations 

Let us finally consider the case £ = £( pt,n t ) where Eq. (14.27P becomes: 


<9 r 7£i 


if /<9£\ 

3(pt + P e ff ) 2 Vdrp/ 


d T p t din t - dip t d T n t 


(4.30) 


The term at the right hand side does not vanish, in general. To lowest order in the gradient 
expansion the diffusion current v a does not contribute to the evolution of n*. Equation (1 4.30j) 
can then be rewritten by using Eqs. (13.21)1 and (13.231) . The result is: 


d T U t 


K 2 N 

3 (Pt + Peff) 



dint dip t 

V dn t J 

[ n t pt + Peff] 


(4.31) 


Similar effects are expected when electromagnetic fields or scalar fields are present together 
with a dissipative fluid (see, for instance, the Erst two papers of Ref. [30] for the case of 
scalar fields and the third paper of Ref. [30] for the case of electromagnetic fields). The 
extension of the present considerations to a multicomponent viscous system is beyond the 
scope of this paper but it is conceptually feasible. 

In summary, the viscous fluids do not necessarily jeopardize the large-scale conservation 
of the curvature inhomogeneities at least as long as £ is a function of the total energy density 
or of the trace of the extrinsic curvature. Conversely the large-scale conservation of the non¬ 
linear curvature perturbations is invalidated whenever £ = £(p t ,n) [or when £ = £(iL, n)]. 
These observations have a counterpart in linear theory which will be discussed in the following 
section. 


5 Back to linear theory 

5.1 Scalar modes in linear theory 

In linear theory variables 7 Z{ and £, defined in Eqs. (14.91) and (14.lip have a well defined limit. 
Let us choose conformally Newtonian frame where the gauge freedom is removed and the 
coordinate system completely fixed 

N 2 (t,x) = a 2 (r)[l + 2<j>(r,x)\, = a 2 (r)[ 1 -2 ^(t, x)\5 ij: (5.1) 

and expand Eqs. (14.9p and (14.lip by assuming that £> and £> are both smaller than one. The 
result of this limit is given by: 


TZi —>■ diTZ, C; —> diC, 


(5.2) 
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implying that 1Z r and Q are, respectively, the spatial gradients of the curvature perturbation 
on comoving orthogonal hypersurfaces and of the density contrast on uniform curvature 
hypersurfaces 


ortnogc 

, i.e@ 


n = -jp- 


uw + 'H<i>) 

U 2 -PL' 


C = -i’ + 


Spt 


3 (Pt + Peff) 


(5.3) 


In this section p t and P e ff denote the background values of the corresponding quantity while 
dp is the first-order fluctuation of the energy density and so on and so forth. In other words 
the conventions will be such that 


p(r,x) = p t (r) + 5p t {r,x), £(r, x) = £(r) + 6£(r, x). (5.4) 


where 5 will denote the first-order fluctuation of the corresponding quantity. The same 
conventions will be employed for all the other variables involved in the discussion. 

The linear order form of the evolution equation for 1Z and ( can be derived in perturbation 
theory and then compared with limit of Eq. (14.211) . The result is is given by the following pair 
of equations where we have included, for the sake of comparison, also the terms that are of 
higher order in the gradient expansion but are consistent with the linearized approximation: 


n' = 

+ 


3P 


a (Pt + Peff) 
3U 2 


Z'VZ + IP) 


U 


Pt + P eff 


dpnad 


SZ + 




a (Pt + Peff) 

C = ^ (C + VO 


-Qt 


3'Hc 2 st 2 

st vV, 


a (pt + Peff) 2£ 2 p(p+ Peff) 

n 


a (Pt + Peff) 

3 n 2 


^Pnad 


K + 0 t 


a (pt + Peff) 


Pt + Peff 

m _i 

- a (Pt + Peff) 3 

£ 


2 t*p(jH + P ef f) 


V 2 (0-^)- 


a (Pt + Peff) 


V 2 ^. 


(5.5) 


(5.6) 


In Eqs. (15.5j) and (I5.6h the term 6 t denotes the three-divergence of the total velocity field. 

Equations (15.5j) and (15.6ft have been derived directly in the linear theory and they repro¬ 
duce the results obtainable by linearizing . Consider then, for the sake of comparison, Eqs. 
(15.5p and (14.271) in the limit 5p na d —> 0. From Eq. (I5.5[) we will have 


n’ 


3H 


a (Pt + Peff) 


m+vo + 


3 V 2 


a (pt + Peff) 


<5£, 


(5.7) 


where 6t has been neglected since it is of higher order in the gradients thanks to the momen¬ 
tum constraint (implying 0 t = SZ 2 (7Z + )/%). 

n The prime will denote a derivation with respect to r and PL = a’/a. 
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To see how tilings work in linear theory we can verify explicitly that the right hand side 
of Eq. (15.7ft vanishes when £ = £(p t ). From the definition of £ in linear theory we have that 
the momentum constraint can be expressed as: 


V 2 0 

C_ + 2 4(p, + P,„Y 


(5.8) 


Neglecting the gradients, Eq. (15.71) is expressible, in the case £ = £(p f ), as: 


K' 


377 rdj\ 

a(Pt + Peff) \dpt) 


pJ(C +VO +'H5p t 


(5.9) 


where we used that <5£ = (d£/dp t )5p t . But now thanks to the definition of £ we have that 
(C + VO = 3 (pt + P e ff)8p t . Thus Eq. (15,9j) implies TV = 0 since, by covariant conservation of 
the background energy-momentum tensor, p[ = — 377(pt + P e ff)- This result does not hold if 
£ = £(tk, n t ) so that, in general, the curvature perturbations induced by relativistic viscous 
fluids are not conserved. This derivation is the linear order counterpart of the discussion 
presented after Eq. (14.271) . 


5.2 Standard gauge-invariant variables 

For infinitesimal coordinate transformations we have that 0 —> 0 = 0 — 77eo — e’ Q and that 
■0 —» 0 = 0 + 77eo- The bulk viscosity coefficient transforms instead as: 

<5£^5£ = 5£-£'e 0 . (5.10) 


Including the gradients in the appropriate entries of the perturbed metric (i.e. Sgij = 
2a 2 (ip5ij — didjE) and Sg M = —a 2 <9*5) and recalling that they transform as B —» 5 = 
B + e 0 — e' and as 17 - > E = E — e the gauge-invariant fluctuations of the bulk viscosity 
fluctuations are: 

~ =5£ + £'(5-5'), (5.11) 

while the Bardeen potentials are, as usual, $ = 0 + (5 — E')' + 77(5 — E') and T = 
"0 — 77(5 — E'). In terms of the explicitly gauge-invariant fluctuations Eq. (15.51) becomes 


= 


377 


a (Pt + ^eff) 

377 2 „ 


+ T) 


77 


Pt + Peff 


8Pnacl 


+ 


£77 


0 / 


377c 


st 


a (p + Peff ) a(pt + Peff ) 2 i 2 p(p + P e ff) 


-v 2 t, 


(5.12) 


where 0* = 9t + V 2 5' is the gauge-invariant expression of the total velocity field. 
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5.3 Fluctuations of the expansion and SJ\f formalism 

We can also compute the total expansion rate with the aim of showing that the presence of 
the bulk viscosity affects the basis of the so-called SAf formalism stipulating that ( can be 
related to the scalar-held perturbations at the initial time computed in the uniform curvature 
gauge once we know the derivatives of the number of efolds with respect to the initial values 
of the unperturbed scalar fields and their derivatives. 

Let us therefore introduce the inhomogeneous generalization of the total number of efolds 

Af(x, t *, Tf) — - j \/ a u a Ndr, (5.13) 

O J T* 

where r* denotes the initial time;all the quantities of the integrand are space-time dependent. 
We recall that in linear theory we can define 

VX = (V/) (0) + 5 (1) (VX) + 5 (2) (VX) + - (5-14) 


Without committing ourselves to a specific gauge choice, from Eq. (15.141) and from the 
fluctuations of the lapse function, Eq. (15.131) becomes: 

Af(x,n,r f ) = JJ(n,T f )+ ^ f f (d t + W 2 E')dr - f f i;'dr, 

O Jr* J t * 

AT(t*,t/-) = f T-idr = f Hdt, (5.15) 

where H = Ha and dt = adr. From the first-order fluctuation of the covariant conservation 
of the total energy-momentum tensor we can express ip' as 

3H 2 


f = \[e t + v 2 £'] + 

3 3 (p t + P e ff) 


l (pt + Peff) 


5£ 


m 


-[6 t + V 2 E'-3(^ + H0)]. 


a (P + Peff ) 

Inserting Eq. (15. 16p into Eq. (I5.15j) we can easily obtain 

3 H 


(5.16) 


JV(x,n,T f ) = Jf(n,r f ) + 1'{ 

£ 


l (Pt + Peff) 


v^n- 


a (Pt + Peff ) 

n 


(.Pt + Peff) 


im+t'w + o] - (c'+vo 


fiPnad r ? 


(5.17) 


where, according to Eq. (15.31) . we used that 5p = 3(£ + ip)(p + P e ff). When 5p na d —* 0 and 
£ = £(/?t) we have: 

7V(f, t *, T f ) = J\T(t*, Tf) (C + ^ dT , (5.18) 
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implying 


A f(x,n,Tf) =JT(n,Tf ) + [£(x,r*) + 0(£,r*)] - [C(f,r/) +'0(f,r / )]. (5.19) 


Let us then evaluate (£ + 0) r/ in the uniform density gauge (i.e. 5p t = 0) while (£ + 0) Tt is 
evaluated in the uniform curvature gauge (i.e. '0 = 0). This choice implies that: 

[C(x, Tf ) + 0(£, Tf)] ->• 0, [£(£, T*) + 0(f, T*)] = 5 +p — y • ( 5 - 20 ) 

Thus, from Eqs. (15.191) and (15.201) we have that 

^ = Af(x, T„, r f ) - A T(n, r f ) = ^ • ( 5 . 21 ) 

6(p t + r'eff) 


The same formulas derived above hold in the single scalar field case (at least up to some point) 
by setting £ = 0 and by recalling that, in the single scalar field case, TZ = —0 — (' H/<p')8<p 
where Sip denotes the scalar field fluctuation. Thus, in the single scalar field case, we will 
have 


[C (f,T*) + 0(f,T*)] = 


TZ* + 0 * + 


V 2 0* 


-» 


—6p. 

p' 


(5.22) 


2 PpiPv+Pv)- 

In analogy with Eq. (15.21)1 . Eq. (15.22)1 can be expressed as 5Af — —(dAf /d<p)*S<p(x,T*), 
which is the standard result of the single scalar field case. 

We can therefore conclude that the presence of bulk viscous stresses affects the 5Af only 
if £ = £(pi, n t ). Conversely, if £ = £(p t ) both TZ and £ are constant and 5AT only depends on 
the values of £ + 0 between the initial and final times characterizing the integrated (total) 
expansion. 


6 Concluding remarks 

The hypothesis of reversibility of the plasma has been dropped by allowing for a relativistic 
dissipative fluid as the dominant source of curvature inhomogeneity in the non-perturbative 
regime. In the general situation the non-linear evolution does not preserve the curvature 
inhomogeneities. However if the bulk viscosity is either a space-time constant or if it de¬ 
pends solely on the total energy density of the plasma, then a further symmetry prevents 
the growth of curvature inhomogeneities. To lowest order in the gradient expansion the 
dependence on the energy density can be traded for a dependence on the trace of the ex¬ 
trinsic curvature. In all these cases the curvature inhomogeneities are non-perturbatively 
conserved. The validity of the linear theory has not been posited as a necessary requirement 
of the derivation. Nonetheless the non-perturbative results have a perturbative counter¬ 
part describable in terms of the conventional gauge-invariant variables namely the curvature 
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perturbations on comoving orthogonal hypersurfaces and the density contrast on uniform 
curvature hypersurfaces. 

The only contribution to the non-perturbative evolution of curvature inhomogeneities 
comes from the bulk viscous stresses while the shear viscosity affects the evolution of the 
traceless part of the extrinsic curvature and it appears to higher order in the gradient ex¬ 
pansion. The only contribution to the inhomogeneous deceleration parameter comes from 
the bulk viscosity. Beyond linear theory the large-scale acceleration can only be driven by 
the bulk pressure and no acceleration can take place thanks to the shear viscosity. 

The results of this investigation confirm that the non-linear conservation of curvature 
perturbations is non a generic phenomenon. However when the curvature inhomogeneities 
are conserved the evolution equations of the extrinsic curvature can be solved in fully inho¬ 
mogeneous terms. These solutions correspond, in the homogeneous and isotropic limit, to a 
quasi-de Sitter stage of expansion. One could therefore speculate that the non-perturbative 
constancy of the curvature inhomogeneities pins down a class of viscous coefficients and, 
ultimately, a specific set of physical properties of the geometry. It is amusing that this logic 
is opposite to the one commonly pursued in linear theory where a quasi-de Sitter stage of 
expansion is postulated to insure, somehow, the perturbative conservation of the linearized 
fluctuations. 
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